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Abstract: We consider a correction to energy levels in a pionic atom induced by the Uehling 
potential, i.e., by a free electron vacuum-polarization loop. The calculation is performed for 
circular states (I — n — 1). The result is obtained in a closed analytic form as a function of 
Za and the pion-to-electron mass ratio. Certain asymptotics of the result are also presented. 
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1. Introduction 

The vacuum polarization effects play an important role in precision theory of energy levels and 
wave functions of atoms. While in conventional atoms such as hydrogen the effects of the self energy 
of an electron dominate over the vacuum polarization, in atoms with a heavier orbiting particle (such 
as a muon or pion) the vacuum polarization is responsible for the dominant QED effects. 

In principle, studies of pionic and other exotic atoms involve effects of the strong interaction. 
There is a direct interaction of the pion and the nucleus. Additionally, the finite-nuclear-size effects 
are important. However, relations between values of various contributions depend on values of the 
nuclear charge Z and of the principal and orbital quantum numbers, n and I, respectively. There are 
certain situations where the vacuum polarization dominates. In a number of other cases, it produces a 
correction, which is not dominant, but still important. 

In this paper we focuss on study of effects of the vacuum polarization within a relativistic frame- 
work. Their leading contribution to the energy levels is related to so-called Uehling potential, a correc- 
tion to the Coulomb potential induced by a free electronic vacuum loop. 

The non-relativistic results in muonic atoms (i.e., the results in the leading order in Za) have been 
known for a while [1] as well as various numerical relativistic calculations for the conventional (see, 
e.g., [2] and references therein) and muonic (see, e.g., [3] and references therein) atoms. 

A relativistic calculation of vacuum polarization effects can be performed analytically. In particular, 
various effects induced by the Uehling potential were studied for the Dirac equation, i.e., for an orbiting 
particle with spin 1/2, in [4, 5, 6]. The results were obtained in a closed analytic form as a function of 
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Za and a mass ratio of the orbiting particle and an electron. However, the result for pionic atoms was 
not derived because it was unclear how to write relativistic perturbative series for particle with zero 
spin. 

Here, we extend our approach onto spinless orbiting particles, such as a pion. Relativistic descrip- 
tion of such a particle is significantly different from that of an electron. A spin- 1/2 particle is described 
by an equation of the form 

where the Hamiltonian iJ is a Hermitian operator and the wave function is its eigenfunction. The 
correction to the energy due to a perturbation by an electrostatic potential SV is of an obvious form 

bEi = {^f\W\^f) . (1) 

In [4, 5, 6] various integrals, which could appear in calculation of such a martix element over the 
Uehling potential in Coulomb bound system were studied. 

The equation, that describes a pion, is Klein-Gordon-Fock equation, which for electrostatic poten- 
tial is of the form 

{{Ei - V{v)f - p^c^ - m^c*] $i(r) = . 

This equation cannot be considered as a problem for eigenfunctions and eigenvalues of a Hermitian 
operator. The solutions of this equation, are not orthogonal and their normaUzation is not trivial. 

While 'exact' solution of the Klein-Gordon-Fock equation with a non-Coulomb potential has been 
possible (see, e.g., study of the problem of falling to the center in [8] and numerical calculations in 
[9]), the perturbative theory with the Coulomb field has not been developed. 

A simple perturbative series for the Klein-Gordon-Fock equation has not been known until now 
and here we apply a recent result obtained in [7]. In particular, the analog of Eq. (1), convenient for 
analytic calculations, is 

This equation is valid for any perturbation that can be described by a certain electrostatic potential and, 
in particular, for the UehUng correction. 

To calculate the Uehling correction we consider the Coulomb problem as an unperturbed case. Its 
eigenvalues and eigenfunctions are known (see, e.g., [10]): 

1 
r 



.s+i / /3r 



nl 

Rniir) = iV(/3r)'^'exp(^-^j iFi(-n + / + l,2s + 2,/Jr) , 
E^^^ = (3) 

where Yim {6, (p) stands for spherical harmonics, iFi (a, b; z) is the confluent hypergeometric function. 
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and m is the mass of the orbiting particle, presumably, a pion. Here, N is a normalization factor, 
unspecified in [10], which is sufficient for our calculation, since the normaUzation is unimportant for 
any application of Eq. 2. 

Our approach can be applied for an arbitrary state. In this paper we consider only circular states 
(l = n — 1), in which the wave function is of the simplest form 

Rniir) = N{l3rye^p(^-^y (4) 

Such states in particular include the Is and 2p states. Some of these states are also of experimental 
interest: the ground state (Is) is the easiest to create. Meantime, high-/ states are easier for theoretical 
interpretation because of a substantial reduction of the strong-interaction effects. 

The Uehhng potential, which is considered as a perturbation of the unperturbed problem with the 
Coulomb potential, is of the form [1 1] 



IT J 1 




_ 2 mec 

where nie is the electron mass. 

The integrals needed for matrix element in (2) are very similar to the ones studied by us previously 
for a Dirac particle [4, 5,6]. The result of calculations of the correction to the energy (2) with potential 
(5) over the wave functions from (3) and (4) is (cf., [4]^) 



SE = (Za) 



2 



(772 + (Za)2) 



2^3/2 



2 



= Jo '^(T^ 



(5 (Za) m 



The results for Kabc{i^) is known in a closed analytic form [4] 



(6) 



KabciK) = B(a + 1/2,1- 6/2 + C/2) (7) 

X 3F2 (c/2, c/2 + 1/2, 1 - 6/2 + c/2; 1/2, a + 3/2 - 6/2 + c/2; k^) 

- ^k''+^ B{a + 1/2,3/2 -b/2 + c/2) 

X 3F2 (c/2 + 1, c/2 + 1/2, 3/2 - 6/2 + c/2; 3/2, a + 2 - 6/2 + c/2; k^) . 

where 3F2 («, 6, c; d,e; z) stands for the generalized hypergeometric function and B(a,h) is the beta 
function. 

The result for a spin-0 particle is similar that a spin- 1/2 particle [4]. We compare the related results 
for the Is state in Fig. 1, where we present results for two values of the mass of the orbiting particle. 



^ Note, that the integral here, Kabc{i^), is slightly different from the one, Iabc{i^, e). used previously in [4, 5, 6]: 
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namely, for m = rUe and m = 250 rrie (for the latter we chose the m value to be somewhere in between 
and m^). The correction is presented in terms of a dimensionless function F{Za) defined as 

The similar plots for 2p and 5g states are presented in Figs. 2 and 3. We note that the results have 
different ranges of convergence, but for physical values Za < 1 the results for excited states are very 
close to each other. 




(a) (b) 

Fig. 1. The Uehling correction to the energy of the Is state for spin-0 and spin-1/2 particles. The parameters are 
(a): m/me = 1 and (b): m/me = 250. 




Fig. 2. The Uehling correction to the energy of the 2p states for spin-0 and spin-1/2 particles. The parameters are 
(a): m/me = 1 and (b): m/me = 250. 

The most important difference is related to the range of values for the nuclear charge Z for which 
the results, derived here for a point-hke nucleus, can be applied. In principle, it is not necessary that 
a correction is well defined for the same range of parameters as the leading term is. For instance, in a 
Dirac atom with a point-like nucleus the hyperfine interval (which can be considered as a perturbative 
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Fig. 3. The Uehling correction to the energy of the hg states for spin-0 and spin-1/2 particles. The parameters are 

(a): m/rn^ — 1 and (b): m/rrie = 250. 



correction to the electron energy levels due to an interaction with the nuclear spin) for the Is state is 
well defined for a narrower range (Za < \/3/2) than the leading term, namely, the Coulomb energy. 
Nevertheless, the Uehling correction for both spin-0 and spin-1/2 particles is well defined in the same 
Z range as the related Coulomb term. In particular, for the Is state of a Klein-Gordon particle the 
condition is Za < 1/2, while for a Dirac particle it is Za < 1. For the hyperfine splitting and the g 
factor of a Dirac particle the situation is similar — the Uehling correction is valid for the same range 
of Za as the Dirac term (for the reference to the Uehling correction in the Is state see [4] for the Lamb 
shift, [5] for the hyperfine splitting and [6] for the g factor of a bound electron). 
The result we obtained above (6) contains two free parameters: Za and 



Za m 

n me 



where 



y/rj^ + {Zaf 
In a pionic atom 



K„ ~ K„ X ( H- {Za) 



2n^{2l + 1) 



while in other mesonic atoms this ratio can be even higher. 

The general expression is rather complicated and we present here some useful asymptotics applying 
an expansion in either of these two parameters. First we consider a case of Za <C 1 which allows a 
non-relativistic expansion. The result is found to be 



n I n[2n — 1) 



+ 



(8) 
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where 



f 

Jo 



dv 



In 



1 + — v'^ 



Various asymptotics of the integral above were studied in [5, 12] (see also below). 

An alternative case with 1 is reaUzed for low-lying states in medium-Z and heavy pionic 

atoms. The result is of the form 



SE 



a [Z(\)-t 



(772 + (Za)2)3/2 



2{Zaf 
2r]-l 



^ln(2;5„) + ^V(l)-^ 



3^^(2^)"3M)^(^^-^)+2Hr^ 

T] r]{27j + l) + 2{Zaf 



+ ... 



(9) 



where tp{x) = ^ ln(r(a:;)) is the logarithmic derivative of the Euler gamma function. 

To check our calculations we have examined a double asymptotic behavior for the case Za <^ 1 
and K„ » 1. The results derived from (8) and (9) agreed. The result for the double asymptotics^ 



5E = - 



a {Za)'^mc^ 



+ 



IT n' 

n{2n - 1) 
3 



2 , X 2/,,^. , x\ 5 Tin n(2n-|-l) 
- ln(2/c„) + - (V(l) - i>{2n)) -- + — - ^ , ^ + , 



1 



2kI 



2n 



— + 1 ln(2Ac„) + -V(l)-- + — 



2n 



- 1 xl}{2n) + 2 V'(2n - 1) - 2n ^'(2n) 



3n 

4;/- + (in - 1 



+ ... 



is more transparent than the single-parameter expansions (8) and (9). 

Other asymptotic cases of interest are related to k„ <C 1. These limit is approached at a pionic 
atom for high n. For arbitrary Za the expansion is 



5E = - 



a {Za)'^mc^ 
^3(r;2 + (Za)2)3/2 



2?H- 1 
2?7 



-(Za)2B(5/2,r?-l/2)K^' 



2r;-l 



+ (r? + l)(r?-2(Za)2)B(5/2,,7)«2'7 

while for small Za we arrive at 

a {ZQ)'^rric^ 



SE = - 
-I- 



TT 3n2 
(2n - 1) 



(n 1)5(5/2, n)4" - n(2n + 3)^(5/2, n + 1/2) k^' 



2n+l 



(10) 



(11) 



(2n-|-l)B(5/2,n- 1/2)k: 



2n-l 



^ Compare the non-relativistic limit (i.e., the leading term in (Za)) to the related result in [4] (corrected accordingly to [6]) 
and [13]. We are grateful to J. Soto for finding a misprint in sign in one of terms in (10). 
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- n(n + 1) (2 \n{Kn) + ^{n) - Hn + 3/2) + ^ 

+ (n + lK2n + 3) )^(^/^'-)-"" + - 

Above we considered the Uehling correction to the energy of the circular states in a pionic atom 
and found the general expression and its various asymptotics. The same approach can be applied for 
any state. The result for low Z, where the non-relativistic effects dominate, is similar to that for a 
bound Dirac particle. At higher Z the correction is different and, in particular, its has a different range 
of applicability (see, e.g.. Figs. 1-3). More detail will be presented elsewhere. 

Acknowledgements 

This work was supported in part by the RFBR under grant 03-02-04029 and by DFG under grant 
GZ 436 RUS 11 3/769/0-1. Part of the work was done during visits by EK to Max-Planck-lnstitut fur 
Quantenoptik and we thank MPQ for their hospitahty. We are grateful to Aleksander Milstein for useful 
discussions. 



References 

1. G. E. Pustovalov, Sov. Phys. JETP 5, 1234 (1957); 

D. D. Ivanenko and G. E. Pustovalov, Adv. Phys. Sci. 61, 1943 (1957). 

2. P. J. Mohr, G. Plunien, and G. Soff, Phys. Rep. 293, 227 (1998). 

3. E. Borie and G. A. Rinker, Rev. Mod. Phys. 54, 67 (1982). 

4. S. G. Karshenboim, Can. J. Phys. 76, 169 (1998); see also: JETP 89, 850 (1999). 

5. S. G. Karshenboim, V. G. Ivanov and V. M. Shabaev, Can. J. Phys. 76, 503 (1998); see also: JETP 90, 
850, (2000). 

6. S. G. Karshenboim, V. G. Ivanov and V. M. Shabaev, Can. J. Phys. 79, 81 (2001); see also: JETP 93, 

477 (2001). 

7. R. N. Lee, A. I. Milstein, and S. G. Karshenboim, Phys. Rev. A 73, 012505 (2006). 

8. V. S. Popov, Sov. Jour. Nucl. Phys. 12, 235 (1971); 
V. S. Popov, Sov. Phys. JETP 32, 526 (1971) ; 

Y. B. Zeldovich, V. S. Popov, Sov. Phys. Uspekhi 14, 673 (1972). 

9. W. Fleischer and G. Soff, Z. Naturforschung 101, 703 (1984). 

10. A. S. Davidov. Quantum Mechanics. Pergamon Press, Oxford. 1976. See Sect. 58. 

11. J. Schwinger. Particles, Sources, and Fields, Vol. 2. Perseus Book, Reading. 1998. 

12. S. G. Karshenboim, U. Jentschura, V. G. Ivanov and G. Soff, Euro. J. Phys. D2, 209 (1998). 

13. D. Eiras and J. Soto. Phys. Lett. B491, 101 (2000). 



2005 NRC Canada 



